Electromagnetic response of a metamaterial with field-gradient-induced transparency 
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We investigate a dynamically controllable electromagnetically induced transparency-like meta- 
material. The unit structure of the metamaterial consists of a low-quality-factor resonator and a 
high-quality-factor resonator. The field gradient of the incident electromagnetic wave in the trans- 
verse direction induces coupling between these two types of resonators and causes a transparency 
phenomenon. We present the simulation and experimental results for the dynamic control of the 
electromagnetic response. 

PACS numbers: 78.67.Pt, 42.25.Bs, 78.20.Ci 



I. INTRODUCTION 

A number of studies have focused on controlling the 
propagation of electromagnetic waves using metamate- 
rials that consist of periodically or randomly arranged 
artificial subwavelength structures. Metamaterials can 
even be used to obtain extraordinary media that do not 
occur in nature, such as media with a negative refractive 
index. 1 Metamaterials have been applied to demonstrate 
various phenomena, such as subwavelength imaging, 2,3 
cloaking^ and no-reflection propagation^ 7 - 

On the other hand, electromagnetically induced trans- 
parency (EIT) has attracted considerable attention in re- 
cent years as a means by which to achieve extremely low 
group velocity propagation of electromagnetic waves; 8,9 
Electromagnetically induced transparency is a quantum 
phenomenon that arises in three-state A-type atoms in- 
teracting with electromagnetic fields, and causes the 
suppression of the absorption of incident electromag- 
netic waves in a narrow frequency range. Based on 
the Kramers-Kronig relations^ steep dispersion, or low 
group velocity, can be achieved in the frequency range. 
In fact, a group velocity of 17m/s ~ 10~ 7 co (where cq is 
the speed of light in a vacuum) has been observed using 
the phenomenon^ 

Since a somewhat complicated arrangement is nec- 
essary in order to achieve EIT, a number of studies 
have focused on mimicking the effect in classical sys- 
tems, such as optical resonators ^r— waveguides^ and 
metamaterials.— ~— However, few studies have focused 
on the dynamic control of the transmission spectrum in 
classical EIT-like systems. In the present study, we inves- 
tigate an EIT-like metamaterial having properties that 
can be controlled dynamically. 

First, we describe the structure of the EIT-like meta- 
material. The unit cell of the metamaterial consists of a 
low-quality-factor (Q) resonator and a high-Q resonator. 
The field gradient of the incident electromagnetic wave in 
the transverse direction induces coupling between these 
two types of resonators and causes a transparency phe- 
nomenon in the metamaterial. Then, we describe how 
to control the electromagnetic response of the metama- 



terial. We present the experimental results for the dy- 
namic control of the transparent frequency using variable 
capacitance diodes and the simulation results for the dy- 
namic control of the frequency width of the transparency 
window. Although the experiment is performed in the 
microwave region, it may be possible to perform similar 
experiments in the terahertz and optical regions. 



II. THEORY 

Figure QJa) shows a capacitively loaded strip 
metamaterial^ The metamaterial can be regarded as 
series inductor-capacitor resonant circuits that consist 
of metal strips in the x direction and gaps between the 
neighboring structures in the x direction. When a plane 
electromagnetic wave is incident on the metamaterial, 
as shown in Fig.QJa), the metamaterial behaves as a 
Lorentz-type dielectric medium because currents are in- 
duced in the x direction by the incident electric field and 
the electric dipoles are excited. 

We now consider the metamaterial shown in Fig.QJb). 
Each capacitively loaded strip is connected to the neigh- 
boring structures in the y direction. The metamaterial 
shown in Fig.[]Jb) includes two types of series inductor- 
capacitor resonant circuits. One consists of a metal strip 
in the x direction and a gap between the neighboring 
structures in the x direction, which is the same as the res- 
onant circuit in the capacitively loaded strips. The other 
consists of a metal strip in the y direction and a gap be- 
tween the neighboring structures in the x direction. (The 
numerically calculated current and electric field distribu- 
tions in the metamaterial are shown in Fig. U and can be 
used to determine the portions of each structure that act 
as inductors and capacitors.) The former resonant cir- 
cuit is referred to as resonator 1, and the latter resonant 
circuit as resonator 2. These two types of resonators are 
magnetically coupled to each other. The x-polarized in- 
cident electric field can directly excite only resonator 1. 
Resonator 2 cannot be excited directly because the di- 
rection of current flow in resonator 2 is perpendicular to 
the direction of the incident electric field. 

Magnetic coupling occurs between the two types of 
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FIG. 1: (a) Metal structure of the capacitively loaded strips 
and (b) the field-gradient-induced-transparency metamate- 
rial. The solid (dashed) arrows represent the current in res- 
onator f (resonator 2). 



resonators when the field gradient of the incident elec- 
tromagnetic fields exists. When the x component of the 
electric field has a gradient in the y direction, the in- 
duced current in resonator 1 becomes y dependent. The 
difference between adjacent currents creates a magnetic 
flux through the loop of resonator 2 and induces anti- 
parallel currents via the electromotive force. When the 
incident electromagnetic wave is a plane wave, this cou- 
pling vanishes due to the uniformity in the y direction of 
the currents in resonator 1. Thus, only resonator I can 
be excited. In this case, the metamaterial behaves as a 
Lorentz-type medium, as in the case of the capacitively 
loaded strips. 

The resonance mode in resonator f is electric dipole 



resonance whereas that in resonator 2 is weak radiative 
resonance, i.e., trapped-mode resonance^ This implies 
that resonator 2 has a higher Q value than resonator 1. 
Therefore, this system is similar to the classical model 
of EIT. 27 When the difference between the frequency of 
the incident electromagnetic wave and the resonance fre- 
quency of resonator 2 is large, resonator 2 is excited only 
slightly and the properties of the metamaterial are de- 
termined primarily by resonator 1. Thus, the metamate- 
rial behaves almost as a Lorentz-type dielectric medium. 
When the frequency of the electromagnetic wave is nearly 
equal to the resonance frequency of resonator 2, resonator 
2 is excited through the coupling with resonator 1. Res- 
onator 2 stores the electromagnetic energy with low loss 
because its Q value is high. The current in resonator 1 
is the sum of the current induced by the incident electro- 
magnetic wave and the current induced by resonator 2. 
Due to destructive interference between these currents, 
the current in resonator 1 might vanish, at which point 
the metamaterial would become transparent. Since the 
transparency phenomenon is induced by the field gradi- 
ent of electromagnetic waves, we refer to the metamate- 
rial as a "field-gradient-induced-transparency metamate- 
rial." 

We can tune three parameters in order to control the 
electromagnetic response of our metamaterial: the res- 
onance frequencies of resonators 1 and 2, and the field 
gradient in the y direction of the incident electromag- 
netic wave. The first two parameters determine the cen- 
ter frequency of the absorption line and that of the trans- 
parency window, respectively. The third parameter cor- 
responds to the coupling strength between resonators 1 
and 2 and determines the frequency width of the trans- 
parency window. 



III. EXPERIMENT ON DYNAMIC CONTROL 
OF TRANSPARENT FREQUENCY 

We fabricated the metamaterial shown in Fig. [2] us- 
ing printed circuit board technology. The thicknesses of 
the copper layer and the FR-4 substrate (relative per- 
mittivity: 4.5) of the printed circuit board were 35 /im 
and 1.6 mm, respectively. We used a structure that was 
slightly different from that shown in Fig.QJb) to adjust 
the resonance frequencies of resonators 1 and 2. We in- 
troduced variable capacitance diodes (Infineon BB857) 
for dynamic control of the resonance frequency of res- 
onator 2. Their capacitances were tuned by applying a 
direct current reverse bias voltage^ 

We measured the transmissivity and the group delay 
of the fabricated metamaterial using a network analyzer. 
A layer of the metamaterial was placed in a rectangular 
waveguide, the cross-sectional dimensions in the x and y 
directions of which were 34.0 mm (height) and 72.1mm 
(width), respectively. The measurement was performed 
in the frequency region in which only the TEio mode 29 
can propagate. The waveguide walls parallel to the yz 
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FIG. 2: Schematic of the metamaterial used for the ex- 
perimental demonstration. The geometrical parameters are 
li — 29.1mm, h = 11.8 mm, I3 — 5.0 mm, Wi = 3.0 mm, 
11)2 = 2.0 mm, W3 = 0.5 mm, IU4 = 1.0 mm, gi = 0.4 mm, and 
g2 — 1.0 mm. 
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FIG. 3: Transmission spectra of the metamaterial shown in 
Fig. [2] with diodes (black) and without diodes (gray). The 
reverse bias voltage on the diodes is zero. 



plane are equivalent to periodic boundaries because the 
electromagnetic fields are uniform in the x direction. The 
width of the waveguide determines the field gradient of 
the electromagnetic waves in the y direction. 

Figure [3] shows the measured transmission spectra for 
the cases in which variable capacitance diodes are con- 
nected (black) and not connected (gray). These two cases 
correspond to the field-gradient-induced-transparency 
metamaterial and the capacitivily loaded strips. In the 
latter case, a simple absorption spectrum is observed, as 
in the case of a Lorentz-type medium. In the former case, 
a transmission window is seen in the absorption band, as 
in the case of an EIT-like medium. Figure @|a) shows 
the numerically calculated current and electric field dis- 
tributions at the resonance frequency of the capacitively 
loaded strips shown in Fig.[IJa). Figure SJb) shows the 
distributions of resonator 2 in the field-gradient-induced- 
transparency metamaterial shown in Fig.QJb). We used 
the COMSOL Multiphysics finite element solver. A sec- 
tion of the metal at the center of the I-shaped structure 
was replaced with a piece of dielectric medium so that we 



can tune the resonance frequency of resonator 1 to that 
of resonator 2. The current distribution and the electric 
field distribution enable us to locate the portions of the 
structure that act as inductors and capacitors. The metal 
strip in the vertical direction (horizontal direction) con- 
stitutes the inductor for resonator 1 (resonator 2), and 
the gap between the neighboring I-shaped structures in 
the vertical direction constitutes the capacitor for both 
resonators. It is confirmed that the electric dipole res- 
onance arises at the resonance frequency in the capaci- 
tively loaded strips and that the trapped-mode resonance 
arises at the resonance frequency of resonator 2 in the 
field-gradient-induced-transparency metamaterial. Note 
that no significant current flows in the vertical direction 
in the case of Fig.@|b). 

Figure [5] shows the measured transmissivity and the 
group delay of our metamaterial as a function of fre- 
quency for three different reverse bias voltages on the 
variable capacitance diodes. The center frequency of 
the transparency window (low group-velocity band) is 
increased by increasing the reverse bias voltage, i.e., by 
decreasing the capacitance of the diodes. 

Next, we present a quantitative analysis of our meta- 
material. The charges q\ and q 2 in resonators 1 and 2, 
respectively, are assumed to satisfy the following equa- 
tions: 

^ + ( 27 r 7 0^ + (27T/0V - (2™) 2 g 2 

= a'Ec- i27Tft + c.c, (1) 

^ + (27r 72 ) ^ + (27r/ 2 ) 2 g 2 - (2™)V =0, (2) 

where / is the frequency of the incident field, /1 (/ 2 ) is 
the resonance frequency of resonator 1 (2), 71 (72) is the 
loss in resonator 1 (2), k is the coupling factor between 
resonators 1 and 2, E is the complex amplitude of the 
incident electric field, a' is a constant, and c.c. stands 
for the complex conjugate of the preceding term.— The 
electric susceptibility of the metamaterial is given as 

m _ -a(f 2 + mf-fi) ( <h\ 

XcU> (f 2 + i-nf - fl)(f 2 + 172/ - /!) -« 4 n/' 

(3) 

where qi = qie l2lx f l + c.c. and a is a constant. As- 
suming that fi / 2 and ji 3> 72, we find that the 
transmittance at / 2 , which is the transparent frequency, 
decreases with decreasing k, and the bandwidth of the 
transparency window can be approximated as k 2 // 2 . The 
group velocity at / 2 is found to decrease with decreasing 
k until k 2 // 2 reaches approximately ^7172 Therefore, 
we can consider the minimum bandwidth of the trans- 
parency window to be ^7172 • The transmissivity t of 
the metamaterial is written as follows^ 

t = ^1 (4) 

[(Z 2 + Zi) 2 e- ifc ^-(Z 2 -Zi)V fc ^]e ifel£i ' 

where k x = (2^//c )^l - [c /(2 W /)] 2 , k 2 = 
(27r//c ) v /l + Xc- [c /(2W)] 2 , Z x = 27r//io/fei, Z 2 = 
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FIG. 4: Current distributions (left panel) and electric field distributions (right panel) in (a) the capacitively loaded strips 
[Fig-[H a )] at the resonance frequency /i of the electric dipole resonance and in (b) the field-gradient-induced-transparency 
metamaterial [Fig.fjjb)] at the resonance frequency /2 of resonator 2. In the right-hand panels, black (white) corresponds to a 
large (small) electric field. The outer boxes represent the waveguide walls (perfect electric conductors). A section of the metal 
at the center of the I-shaped structure is replaced by a piece of dielectric medium in order to tune /i and /2 . The relative field 
strengths in (a) and (b) are not scaled. 



27r//^o/A:2r^ w is the waveguide width, fj,Q is the per- 
meability of the vacuum, and d is the thickness of the 
metamaterial. We adopt d = 1.7 mm from the physical 
thickness of the metamaterial. It is difficult to define the 
thickness of a single-layer metamaterial. However, in the 
present case, the change in d has little effect on the fol- 
lowing fitting results. We compare t with the measured 
transmission spectra. 

The fitted curves of and the fitting parameters are 
shown in the upper row of Fig. [5] Here, we adopt the val- 
ues of /i = 2.93 GHz and 71 = 0.20 GHz derived from the 
measured transmissivity and group delay (not shown) of 
our metamaterial without the diodes and assume that /2, 
72, k, and a are fitting parameters. Although the mea- 
sured and fitted transmissivities are slightly different in 
the lower-frequency region, they agree rather well with 
each other. The group delay {d[arg (i)]/d/}/27r calcu- 
lated using the obtained fitting parameters is shown in 
the lower row of Fig.[SJ The measured group delay is in 
good agreement with the calculated value. 

Based on the above results, the condition whereby 
K 2 /f2 > 7i > 72 is satisfied. This implies that low-Q 
resonator 1 is strongly coupled to high-Q resonator 2. In 
addition, f 2/^/2 is found to increase with the reverse bias 
voltage. This is because the resistance of the variable ca- 
pacitance diode, which determines the loss of resonator 2, 
decreases with increasing reverse bias voltage. The value 
/2/K is independent of the reverse bias voltage, i.e., the 



coupling strength between resonators 1 and 2 does not 
depend on the resonance frequency of resonator 2. A 
possible reason for the slight difference between the mea- 
sured and fitted values is that we used one layer of the 
metamaterial in the experiment whereas in the theoret- 
ical model, we assumed a homogeneous medium. The 
oblique propagation in the waveguide^ may also have to 
be considered. 

In order to confirm the Q values of the two types of res- 
onators, we performed COMSOL simulations and found 
that Qi = /1/71 = 9.1 and Q 2 = / 2 /7 2 = 71. The dielec- 
tric loss tangent of the substrate was assumed to be 0.02. 
The Q values calculated from the fitting result for the 
reverse bias voltage of 30 V were Q\ = 15 and Q2 = 83. 
The experimental values are slightly larger than the sim- 
ulated values for both types of resonators. This might 
be because an actual value of the dielectric loss tangent 
of the substrate was smaller than 0.02. The ratio Q2/Q1 
for the experiment might be smaller than that for the 
simulation because of the resistance of the diodes. 



IV. SIMULATION OF DYNAMIC CONTROL OF 
COUPLING STRENGTH 

We performed numerical simulations using COMSOL 
Multiphysics to confirm that the coupling strength be- 
tween resonators 1 and 2 depends on the field gradient 
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FIG. 5: Transmissivity (upper row) and group delay (lower row) of the fabricated metamaterial as a function of frequency. 
Black crosses and gray solid curves represent the measured data and the fitting curves, respectively. The reverse bias voltage 
on the variable capacitance diodes are (a) OV, (b) 10.5 V, and (c) 30 V. 
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FIG. 6: Dependence of the transmission spectrum on the 
waveguide width or the field gradient in the y direction of 
the incident electromagnetic wave. The waveguide width for 
the dashed curve is 2.2 times larger than that for the solid 
curve. 



Wb = 4.0 mm, and w c — 0.5 mm, and the thickness of 
the metal was 0.5 mm. The metal pattern was on a di- 
electric substrate of 0.5 mm in thickness having a rela- 
tive permittivity of 4.5. In order to tune the resonance 
frequency of resonator 1 to that of resonator 2, we re- 
placed part of the metal at the center of the I-shaped 
structure with a piece of dielectric medium, the length 
and relative permittivity of which were 1.8 mm and 24.5, 
respectively (see Fig. [4]). The dimensions of the cross sec- 
tion of the waveguide in the x and y directions when the 
solid curve (dashed curve) was obtained were 34.0 mm 
and 70.0 mm (34.0 mm and 153.7 mm), respectively. We 
confirmed that the transparency window narrows as the 
waveguide width increases. In other words, the smaller 
the field gradient in the y direction of the incident electro- 
magnetic wave, the weaker the coupling strength between 
resonators 1 and 2. 



in the y direction. We calculated the transmissivity of 
the metamaterial placed in a rectangular waveguide. We 
changed the width of the waveguide in order to control 
the field gradient in the y direction of the incident elec- 
tromagnetic wave and calculated the transmissivity for 
several values of the waveguide width. 

Figure [6] shows the transmission spectra of the meta- 
material shown in Fig.[ljb) for the TEio mode. The 
geometric parameters of the metamaterial were Z a = 
25.1mm, 1^ = 18.0 mm, l c = 2.4 mm, w a = 4.0 mm, 



V. SUMMARY AND DISCUSSION 

We investigated a dynamically controllable metama- 
terial that behaves as an EIT-like medium. The unit 
structure of the metamaterial is composed of a low-Q res- 
onator (resonator 1) and a high-Q resonator (resonator 
2). When a plane electromagnetic wave is incident on the 
metamaterial, only resonator 1 can be excited. Resonator 
1 couples to resonator 2 when the incident electromag- 
netic wave has a field-gradient in the transverse direction, 
such as a Gaussian beam. The system is similar to the 
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classical analog of EIT, and the metamaterial behaves as 
an EIT-likc medium. We described a method by which 
to control the electromagnetic response of the metamate- 
rial. We experimentally demonstrated dynamic control of 
the transparent frequency, i.e., the resonance frequency 
of resonator 2 of the metamaterial, using variable ca- 
pacitance diodes. The simulations confirmed that the 
frequency width of the transparency window, or the cou- 
pling strength between resonators 1 and 2, can be con- 
trolled by the field gradient in the transverse direction of 
the incident electromagnetic wave. 

In our experiment, the ratios of the wavelength to the 
length of the unit structure in the x and y directions are 
2.3 and 6.5, respectively. These ratios are larger than 2, 
and thus the diffraction can safely be neglected. For nor- 
mal incidence, in particular, the diffraction does not arise 
when these ratios are larger than 1. This implies that 
effective-medium parameters of the metamaterial can be 
defined. That is, the metamaterial can be regarded as a 
homogeneous medium. For a larger ratio in the x direc- 
tion, we could reduce the length of the unit structure in 
the x direction while preserving the resonance frequencies 
of resonators 1 and 2 by reducing and w\. 

Although we performed our experiment in the mi- 
crowave region, it is not difficult to conduct similar ex- 
periments in higher-frequency regions, such as the ter- 
ahertz and optical regions. Since the structure of our 
metamaterial is planar, we can fabricate the metamate- 
rial in the optical region by standard photolithography 
techniques. We measured the transmission spectrum of 
the fabricated metamaterial in the waveguide, where the 
beamwidth of the electromagnetic wave is similar to the 
wavelength. In the optical region, it is possible to obtain 
a focus spot size close to the diffraction limit, and hence, 
it is possible to control the coupling strength between 
resonators 1 and 2 in free space. We controlled the reso- 
nance frequency of resonator 2 by using variable capaci- 
tance diodes. The resonance frequency can be controlled 
by using the photoexcitation of carriers in a semiconduc- 
tor in the terahertz region^ and the phase transitions 
of liquid crystal 3 ^ and vanadium oxide 33 in the optical 
region. 

Field-gradient-induced-transparency metamaterials 
have numerous applications. The difference between the 



capacitively loaded strips shown in Fig.[lja) and our 
metamaterial shown in Fig.[TJb) is whether or not the 
I-shaped metal structures are connected to neighboring 
structures. Therefore, we can dynamically switch from 
a Lorentz-type medium to an EIT-like medium and vice 
versa by controlling the resistance of the connection 
using the photoexcitation of carriers^! or Schottky 
junction^ This indicates that a dynamic field-gradient- 
induced-transparency metamaterial could be used as an 
amplitude modulator. Since the transparency window 
is narrow and the transparent frequency depends on 
the inductance and capacitance that form resonator 
2, our metamaterial can be used for high-sensitivity 
sensing. The metamaterial can also be applied to 
transverse-mode filters and measurements of beamwidth 
because the transmission spectrum depends on the 
field gradient in the transverse direction of the incident 
electromagnetic wave. 

Field-gradient-induced-transparency bears a resem- 
blance to nuclear quadrupole resonance. In nuclear 
quadrupole resonance, the energy level of nuclear spin 
splits due to the interaction between the electric 
quadrupole moment of atomic nuclei having a nuclear 
spin greater than 1/2 and the electric field gradient at 
the atomic nuclei that is caused by electrons. The split- 
ting width of the energy level increases with the electric 
field gradient. In field-gradient-induced transparency, 
the frequency separation between the two dips in the 
transmission spectrum increases with the field gradient in 
the transverse direction of the incident electromagnetic 
waves. 
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